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Abstract. We propose a conjectural formula for correlation functions of the 
Z-invariant (inhomogeneous) eight-vertex model. We refer to this conjecture as 
Ansatz. It states that correlation functions are linear combinations of products of 
three transcendental functions, with theta functions and derivatives as coefficients. 
I The transcendental functions are essentially logarithmic derivatives of the parti- 

tion function per site. The coefhcients are given in terms of a linear functional TrA 
on the Sklyanin algebra, which interpolates the usual trace on finite dimensional 
representations. We establish the existence of Tr^ and discuss the connection to 
the geometry of the classical limit. We also conjecture that the Ansatz satisfies 
the reduced qKZ equation. As a non-trivial example of the Ansatz, we present a 
new formula for the next-nearest neighbor correlation functions. 
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O . 1. Introduction 

O I Exact description of correlation functions and their analysis is one of the central 

^ I problems of integrable lattice models. Significant progress has been made over the 

last decade toward this goal. In the study of correlation functions, a basic role is 
played by a multiple integral representation, first found for the archetypical example 
Qh! of the spin 1/2- XXZ chain [T71 UHl 120] • Subsequently it has been generalized in 

^ I several directions, to incorporate an external field [201; unequal time non-zero 

temperature ^3] and finite chains ^Hl- Earlier in the literature, extension to elliptic 
models has also been pursued. The free field construction used in the XXZ model was 
rS I extended in [21] to the SOS models, resulting in an integral formula for correlation 

d \ functions of the ABF model. In [23] an integral formula was obtained for the eight- 

vertex model by mapping the problem to the SOS counterpart. A novel free field 
representation of the eight- vertex model is being developed in [23 131] ■ 

Recent studies have revealed another aspect of these integrals. Through examples 
at short distance, it has been observed in the case of the homogeneous XXX chain 
that the relevant integrals can be evaluated in terms of the Riemann zeta function at 
odd integers with rational coefficients [3] . Similar calculations have been performed 
for the XXZ chain |23t i37j . This phenomenon was explained later through a duality 
between the qKZ equations of level and level — 4 [UljT] . Motivated by these works, 
we have established in our previous papers [S] H] an algebraic representation (in the 
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sense no integrals are involved) for general correlation functions of the inhomoge- 
neous six-vertex model and its degeneration ^. The aim of the present paper is to 
continue our study and examine the eight-vertex model. 

We formulate a conjectural formula for correlation functions (the Ansatz) along 
the same line with the six-vertex case. Consider the eight-vertex model where each 
column i (resp. row j) carries an independent spectral parameter ti (resp. 0). The 
object of our interest is the matrix 

hn (tl ; ■ ■ ■ ytn) 
1 3 

ai,- - ,a„=0 

where (■ ■ ■ ) denotes the ground state average in the thermodynamic limit, o"° = 1, a" 
(1 < a < 3) are the Pauli matrices, and T stands for the matrix transpose. Regard 
hr, as a vector via the identification End((C2)®") ~ (C^)®^™^ ^^^^ denote the 

vector corresponding to the identity. Our Ansatz is that hn can be represented in 
the form 

- = 2-"exp Sn- 

\i<j a=l J 

Here uja{t) axe scalar functions given explicitly in terms of the partition function per 
site (see (I2.32p below). The matrices XaJ^i^ are expressible by theta functions and 
derivatives. Leaving the details to Section let us comment on the latter. 

In the six-vertex case, XaJ^^ are defined in terms of a 'trace' of a monodromy 
matrix. Here 'trace' means the unique linear functional 

TrA : f/,(s[2) ^C[g±"]©AC[g±"], 

which for A G Z>o reduces to the usual trace on the A-dimensional irreducible 
representation of t/q(sl2). In the eight- vertex case, we need an analogous functional 
Ttx, defined on the Sklyanin algebra and taking values in the space of entire functions 
involving A, theta functions and derivatives. Compared with the trigonometric case, 
the existence of Tr^ is more difficult to establish. We do that by considering the 
classical limit and showing that, for generic values of the structure constants, the 
computation of the trace of an arbitrary monomial can be reduced to that of seven 
basic monomials. We have also compared our formula for Tr^ with the results by 
K. Fabricius and B. McCoy |Sj for A = 3, 4, 5. In the classical limit, the Sklyanin 
algebra becomes the algebra of regular functions on an algebraic surface in C^, which 
turns out to be a smoothing of a simple-elliptic singularity of K. Saito [23 • The 
reduction of the trace is closely connected with the de Rham cohomology of this 
surface (see Appendix E)). Although we do not use Saito's results for our immediate 
purposes, we find this connection intriguing. 

^Correlation functions of the XXZ and XXX chains are given in the hmit where all the inhomo- 
geneity parameters are chosen to be the same. However we have not succeeded in performing this 
homogeneous limit. 
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In the trigonometric case, it was shown [4j that the functions given by the above 
Ansatz satisfies the reduced qKZ equation. The steps of the proof carry over straight- 
forwardly to the elhptic case, except for one property (the Cancellation identity). 
Unfortunately we have not succeeded in proving this last relation. It remains an 
open question to show that our Ansatz in the elliptic case satisfies the reduced qKZ 
equation. 

To check the validity of the Ansatz, we examine the simplest case n = 2. In this 
case an exact answer for the homogeneous chain is obtained as derivatives of the 
ground state energy of the spin-chain Hamiltonian. Our formula matches with it. 
We also present an explicit formula for the correlators with n = 3. It agrees well 
numerically with the known integral formulas of [221 i2n| • 

The plan of the paper is as follows. In Section 2, we introduce our notation and 
formulate the Ansatz for correlation functions. In Section 3, we discuss the validity 
of the Cancellation Identity and give arguments in its favor. Section 4 is devoted to 
the examples for correlators of the nearest and the next nearest neighbor spins. We 
also discuss briefiy the trigonometric limit. In Appendix El we prove the existence 
of the trace functional. As was mentioned above, the classical limit of the Sklyanin 
algebra is related to an affine algebraic surface, and the the trace functional tends 
to an integral over a certain cycle on it. We explain the connection between this 
picture and K. Saito's theory. In Appendix ^ we give an explicit description of the 
integration cycles. Appendix O contains technical Lemmas about the trace. Finally 
in Appendix O we discuss the transformation properties of the matrices Xa,n\ 



In this section we introduce our notation and formulate the Ansatz for correlation 
functions of an inhomogeneous eight- vertex model, following the scheme developed 
in m. 

2.1. R matrix. We consider an elliptic R matrix depending on three complex pa- 
rameters t,ri,r. We assume Imr > and rj ^ Q + Qt ^. We will normally regard 
r],T as fixed constants and suppress them from the notation. Let 9ait) = 6'a(t|r) 
(0 < a < 4, ^4(t) = 6^0 (t)) denote the Jacobi elliptic theta functions associated with 
the lattice Z + Zr [Tij. We set 



The R matrix is given by 



2. Ansatz for correlation functions 



(2.1) 



R{t): 



Pit) 



e End{V®V), 



[t + v] 



(2.2) 




where V = Cf + © Cf _ . 



'Later on we also assume that 77 is generic 
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The matrix r{t) = Ti2(t) is the unique entire function satisfying 
ri2(0) = Pi2, 

(Tta^ruit) = ruit)a'lcT^ (a = 1, 2, 3), 

ri2(t + ^) = -alru{t)al, 

r,^(^t + ^)= -crfruitWl x e-2-(2t+,+r/2)_ 

Here P G End(V ® V") signifies the transposition Pu ^ v = v ^ u. As is customary, 
the suffix of a matrix indicates the tensor component on which it acts non-trivially, 
e.g. o-f = a" ® 1, = 1 ® a". 

The normahzing factor p(t) is chosen to ensure that the partition function per site 
of the corresponding eight- vertex model equals to 1. Its explicit formula depends 
on the regime under consideration, and will be given later in (|2.28|) . (|2.29p . In each 
case it satisfies 

p{t)p{-t) = l, p{t)p{t-r,)- 



We will often write tij = ti — tj. The basic properties of R{t) are the Yang-Baxter 
equation 

(2.3) -Rl2(i^l2)-Rl3(il3)-R23(^23) = -^23 (^23)-Rl3 (i^l3)-Rl2 (^12) , 

and 

(2.4) R{t) = PR{t)P, 

(2.5) Ri-v) = -2V^, 

(2.6) R,2{t)R2ii-t) = l, 

(2.7) Pi2(t)P2"3 = -Rrsi-t - v)V2s. 

In (12. 5|) . = (1 — P)/2 denotes the projection onto the one- dimensional subspace 
spanned by 



We will use also 



s := f _i_ ® f- — ® f + G V ® V. 



R{t) = PR{t). 



2.2. Sklyanin algebra. Along with the R matrix, we will need the L operator 
whose entries are generators of the Sklyanin algebra jHSl EE] ■ 

Recall that the Sklyanin algebra A is an associative unital C-algebra defined 
through four generators Sa (a = 0, 1, 2, 3) and quadratic relations 

(2.8) [5*0, Sa] = iJhc{SbSc + ScSh)-, 

(2.9) [S,,S,]=i{SoSa + SaSo), 
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where (a, 6, c) runs over cyclic permutations of (1,2,3). The J^c are the structure 
constants given by 

(2.10) oM^e^^M^ 



(2.11) Ja 



where 

£2 = -1, ea = l (a^ 2). 

Here and after, theta functions without arguments stand for the theta zero values, 
e„ = ^„(0)and 9[ = 9'M. 

Since the defining relations are homogeneous, ^ is a Z>o-graded algebra, A = 
®n>o-A.n, where the generators Sa all belong to Ai. We have also a Z2 x Z2-grading, 
A = ©(m,n)eZ2xZ2-^*-™''"\ defined by the assignment Sa E A°', where 

(2.12) = (0,0), 1 = (1,0), 2 = (1,1), 3 = (0,1) GZ2XZ2. 

To make distinction, the Z>o-grading and the Z2 x Z2-grading will be referred to as 
'degree' and 'color', respectively. Thus Sa has degree 1 and color a. 
There are two central elements of degree 2 and color 0, 

3 3 

(2.13) Ko:=Y,Sl K2:=Y,JaSl. 

a=0 a=l 

We call them Casimir elements. 

Introduce the generating function (L operator) 

L{t) ■■=\Y. \'^'^^\'^^ Sa®a- G^®End(\/). 

The defining relations ()2.8|1 . ()2.9|1 are equivalent to 

(2.14) Ruit - s)L,{t)L2is) = L2is)L,{t)Ruit - s). 
From ()2.13|) we have 

(2.15) L, (-jL2(- - v)V,2 = -4 [ ^2 ^0 + ^,K2 j V,2. 

We will be concerned with representations in series (a) of , which are analogs 
of finite-dimensional irreducible representations of s[2- For each non- negative integer 
k, let V'-'^'' denote the vector space of entire functions f{u) with the properties 

fiu + 1) = fiu) = fi-u), f{u + r) = e-^-^^(^"+^)/(«). 
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We have dim V^'^-' = k + 1. The following formula defines a representation n^''^ : A 
End(V('=)) |36j: 3 

(2.16) {n^'\S^)f) (u) 



u) 



Here y/e^ = i, y/e^ = 1 {a ^ 2), and (e=^''^"/)(M) = f{u±i]). In particular, if /c = 1, 
then in an appropriate basis we have 



TT 



=a", (7r(i)®id)L(t) =r(t). 



On V*-'"'', the Casimir elements Kq, K2 act as scalars Ko^k + l), K2{k + 1) respectively, 
where 

(2.17) /.„(A)^4^. AVA)^/'(^'' + ")*'(^"-" 



«i(2.;)' " ' 9i(2.,) 



2 



2.3. The functional Tr^. In order to formulate the Ansatz, we need to consider 
the trace try{fc)7r*^'^')(A) of an element A G ^ as a function of the dimension k + 1. 
The precise meaning is as follows. 

For each A & A one can assign a unique entire function TtxA in A with the 
following properties: 

(i) TrAA|^^^^^= trv(fc)7r('=)(A) holds for all k e Z>o, 

(ii) If y4 G An, TtxA has the functional form 



(2.18) TrtA = 0i(t)"x 



QAflit) ( n: odd ), 

QaA^) - ^9A,2it) ( n: even ), 



where gAfi{t), gA,2(t) and gA,3(t) ■= 9A,i(t + T) - gA,i(t) are elliptic functions 
with periods 1, r. In addition, gA,i{t + 1) = gA,iit)- 
For example, 

(2.19) TrAl = A, 



(2.20) TrA^„ = 26, 



(2.21) TrA^^ = ^Tq^ (FciiXv) ' AF„2(Ar^)) , 
where 

(2.22) F„i(t) = ej^+iir]^^ (^,+i(t + v)Oa+i{t - v)) , 

ot 

(2.23) F„2(t) = ea0a+l{V?-^ {da+l{t + Tl)9o.+l{t - T])) 

or] 

For reference we set 

(2.24) F„3(t) = ej^+i{r]fe^+^{t + r])eo,+i{t - r]). 



^We have modified eq.(6) of [321 by a factor 9i{2ri). 
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Tr^ satisfies also 

(2.25) Ttx{AB) = Ttx{BA), 

(2.26) Tlx (KiA) = K,iX)TrxiA) = 0, 2), 

(2.27) TixA = (AG^(™'"\(m,n)^(0,0)). 

The derivation of as well as (IT^ . ^TT^ is sketched in Appendix O In Ap- 

pendixIXJwe show that, for generic r], any element A G A/ Yla=o\.^o:^ ^] '^^^ writ- 
ten as a €[^0, -^"2] -linear combination of seven monomials: 1, Sq, Si, S2, 83, Sq, S^. 
Hence TtxA is completely determined by the property (j2.25j) - (|2.27|) along with 
(|2.19p - (j2.2ip . However an effective algorithm for the reduction is not known to 
us. The situation is in sharp contrast to the trigonometric case, where a simple 
recursive procedure for calculating Tr^ is available (see |1]). It would be useful if 
one can find a more direct expression for the trace using Q-operators, as is done for 
the XXZ model in 

In this connection, notice that by the above rule, TtxA has a simpler structure 
for elements of odd degree than those of even degree, since in the former case it is a 
polynomial of 61 {t + c) (c G C) without involving derivatives. The subtle difference 
between finite spin chains with odd length and those with even length has been 
noticed in the context of Q-operators [771 1^ 

2.4. The Ansatz. Consider an inhomogeneous eight vertex model, where each col- 
umn i (resp. row j) carries a spectral parameter ti (resp. 0). The Boltzmann weights 

are given by the entries Re\\el{ti) of the R matrix ()2.H) . We choose the normalizing 
factor pit) in ()2.1|) in accordance with the two regimes 

(i) rj,t E iM, —irj > (disordered regime) 

(ii) r], t G M, < (ordered regime) 

In the disordered regime, p{t) = p'^**(t) is given by 

(2.28) p-(t):=e--x^(2,-2t)7(4, + 2t) 



7(2ry + 2t) 7(4r/ - 2t) ' 
7(m) = T{u,Ar],T), 



where 



00 



1 _ g27ri{(j + l)r+{fc+l)(T-«) 

r(M,a,r) := [[ 



^ g27ri(jT+fc(T+M) 

j,k=0 



is the elliptic gamma function j2H]. In the ordered regime, the formula for p(t) = 
p°^'^{t) is changed to 

(2.29) p"'''^{t) = e-^"^''*/>'(t'), 

where p'{t') is given by the right hand side of ()2.28|) with t, 77, r being replaced by 

(2.30) t' = -, v'=-, r' = --, 

T T T 

respectively. 
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By correlation functions we mean the ground state averages (a"^ • ■ -cr"") of a 
product of spin operators on consecutive columns 1, . . . ,n on a same row of the 
lattice. The thermodynamic limit is assumed. We arrange them into a matrix 

hn (tl ) ■ ■ ■ ytn) 
1 ^ 

ai,--- ,an=0 

where (cif )^ = eaCif stands for the transposed matrix. Because of the 'Z-invariance' 
P], it does not depend on ti with i < 1 ot i > n. When ti = ■ ■ ■ = tn, each 
(o""^ ■ ■ ■ 0""") is a correlation function of the infinite XYZ spin chain 

oo 
j=-oo 

at zero temperature (the coefficients will be given below in (j4.5p ). The hn niay 
be viewed as the density matrix of a finite sub-system of length n, regarding the 
rest of the spins as an environment. 

From now on we fix n, and write j = 2n— j + We regard hn as a 2^"-dimensional 
vector through the identification End(V"®") ~ V^^^ given by 



n 

(2.31) Ee„e, ® ■ ■ ■ ® ^ (J[ ej) Ve, ^ ■ ■ ■ (S) Ve„ (g) V^g^ ®---®V^ 

i=i 

where Ee,e> = (4,a4',/3)«,^=±- 

Let us explain the constituents which enter the Ansatz. 

First, we define three functions in terms of the factor p{t) (given in ()2.28p or 
(Eli) by 

d d d 

(2.32) ijiit) :=-log<f{t), uj^it) := g^^ogif{t), Us{t) := —logif{t), 

where we have set 

They are a meromorphic solution of the system of difference equations 

uJi{t-r]) +uJi{t) = qi{t), 

^2{t -V)+ ^2{t) -uJi{t-T]) = q2{t), 

where 

d d d 

(2.33) gi(t):=-logV'(t), g2(t) := ^ log^(t), q3{t) := —logtP{t), 

and 
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The next ingredient are the matrices Xajn^ (a = 1, 2, 3, I < i ^ j < n). Consider 
a 'transfer matrix' 

(2.34) X„(ti,--- 

1 rr. 



We used the functional TrA introduced in the previous section, and 

-^i, ^ (^) ^1 ) • • • ) ^n) 

:= L2{t -t2 -!])■■■ Ln{t -tn- r])Ln(t - t„) ■ ■ ■ L2(t - 12). 

Notice the presence of the permutation P12 and the projectors "^^^'^22 ()2-34p . 
For i < j, we define 

(2.35) := Mj^'"'''(ti, ■ ■ ■ , tn)Xn{ti, tj, ti, ■ ■ ■ ,ti, ■ ■ ■ ,tj, ■ ■ ■ , t^) 





X I 


R(:'^-)(ti,--- ,t„)-\ 




Here K 


in''^^ stands for the product of R matrices 




(2.36) 


^'■^■^(ti,--- ,tn: 


) 






'■= Ri,i-l{ii,i-l) 


■ ■ ■ -R2,l(^i,l) 






^ Rj,j-^ih,j-^) 


■ ■ ■ -Rj+2,i+l(ij,i+l) ■ 


Ri+l,i{tj,i~l) ■ ■ ■ ■ -R3,2(^j,l) 




X Rjrj-iiti-i^) 


■ ■ ■ Rl2{'tl,i) 






X Rj^-jitj^ij) 


' ■ ■ ■ -^m,i+2(^i+i,j) ■ 


■ ■ ■ R2,3{tl,j)- 


Finally, for all i ^ j, Xa]n are defined by 




(2.37) 




:=X(^'^-)(ti,--- ,t„) 


-U,A^^Xi^'^\h,--- ,tn), 


(2.38) 


'^^2,n '*(^l5 ■ ■ ■ ; ^n) ■ 


= -r/A«X^-)(ti,. 


; ^n) ? 


(2.39) 


'^^3,n\tl^ ■ ■ ■ ,tn) 


:= A«X(*'^■)(tl,■■■ 


,t„)-rAPx(^'^')(ti,--- ,t„ 


where 


c = ~29[/9i{2r]) and 








A«/(---,t„--- 


) = /(■ ■ ■ + a, ■ ■ 


■)-/(■■■, i^,---)- 


As we 


show in Appendix IHl the Xa,n{ti, ■ ■ ■ , 


t„) are doubly periodic in 



periods l,r. The only exception is the case a = 1, k = i or j and with respect to 
the shift by r, where the transformation law becomes 

A^^gn^i, ■ ■ ■ ,tn) = ^^l^ifih, ■ ■ ■ ,t„) = Xi'fih, . . . ,Q. 
Conversely we have 

xi^'Kh, ■ ■ ■ , t„) = c (xi:^ (ti, ■ ■ ■ , t J - ^x(:;^'^(ti, ■ ■ ■ , 
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We compute the trace in the formula ()2.35|) by using the formulas ()2.17|) and ()2.2H) . 
The separation of Xn'-'^ into two parts xf^-* and Xg*^"* comes from that of Tr^S^ 

into F„i(Ar/) and F„2(Ar/). Note that xj^;'^ = x['f and X^,i^ = -X^^ for a = 2, 3. 
We are now in a position to state our conjecture. Let 

n 
p=l 

be the vector corresponding to the identity by the map (j2.3ip . 

Conjecture. Correlation functions of the inhomogeneous eight- vertex model are 
given by the formula 

(2.40) hnih, ■■■ ,tn)= 2-"exp (j2J2uJa{Uj)Xi]'^\ti, ■ ■ ■ ,t„) j s„, 

\ i<j a=l / 

where Ua{t) and X^a,n are defined respectively by (ITT?^ and ^IMi - ^I^ . □ 

3. Reduced qKZ equation 
The hn is known to satisfy the following set of equations [T^ : 

(3.1) hn{- ■ ■ , tj^i,tj, ■ ■ ■ ) 

= Rjj+i{tjj+i)Rj^j(tj+ij)hn{- ■ ■ , tj,tj+i, ■■■) (1 < J < n - 1), 

(3.2) hn{---,tj-r],---) = Al^\h,---,tn)hn{---,tj,---), 

(3.3) V^j ■ hniti, ■ ■ ■ , tn)l,...,n,n,...,l = -^Silhn-l{t2, ■ ■ ■ , tn)2,...,n,n,...,2- 

Here 

(3.4) A(f)(ti,-- - ,y 

In this section, assuming a conjectural identity, we explain that these relations are 
valid also for the Ansatz. 

3.1. Properties of ^In'''^ . Consider the expression 

3 
a=l 

which enters the Ansatz ()2.40|) . In for the XXZ model, the following relations 
are derived. 

Exchange relation: 

(3.5) Rk,k+l{tk,k+l)Rk+T;ki^k+l,k)^n'''\' ' ' ,tk,tk+l, ■ ■ ■) 

= fi(7'=W'''*(^))(- ■ ■ • ■ ■)Rk,k+i{tk,k+i)Rk+ij^{tk+i,k), 

Here vr^ signifies the transposition (A;, k + 1). 
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Difference equations: 

(3.6) ,tk-Vr-- ,tn) 

= AW(ti, ■ ■ ■ ,Q&:^^\t,, ■ ■ ■ ,QAl'\t,, ■ ■ ■ ,t„)-i (k ^ 

= A«(ti, ■ ■ ■ ,t„) ■■■,tn)+ Yj^'''\h, ■ ■ ■ ,t„)) Sn. 

In the last line, we have set 

3 

a=l 

where qa(t) are given by (j2.33p . 
Recurrence relation: 

(3.7) V--,Q^:'^\h,--- ,Q 
{1 = i < j < n), 

Commutativity: For distinct indices i,j,k,l, 

(3.8) fi(*'^)(ti,--- ,t„)fi^')(ti,--- ,t„) = ^('=''^^1^ 
Nilpotency: 

(3.9) ,tn)nl'^^\h,--- ,tn) = o if Kj}n{fc,/}^0. 

The proof of these relations given in are based only on the properties ()2.3|) . ()2.4p - 
(j2.7|) of the R matrix and (j2.14j) . (j2.15j) of the L operator. Hence they carry over to 
the elliptic case as well. 

As is shown in |4J, Proposition 4.1, the equations (j3.5|) - (j3.9|) guarantee the va- 
lidity of the fundamental properties ()3.1|) . ()3.2|) . ()3.3|) for the Ansatz, provided one 
additional identity holds: 

Cancellation identity: 

(3.10) ■■■,tn)+ (A«(ti, ■ ■ ■ ,t„)-l - 1) j S„ = 0. 

So far we have not been able to prove the cancellation identity. In the next 
subsection, we suggest a possible approach toward its proof. 

3.2. Cancellation identity. Set 

(3.11) Q«(ti,---,tO 

■ ■ ■ ,tn) + K'Hii, ■ ■ ■ ,t^r' - 1) ^n. 

We regard it as a matrix via the isomorphism ()2.3ip . 
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Besides the obvious translation invariance, Qn = Qn^ has the following properties. 

(3.12) n"=2 ^i(2ii,i) ■ Qniti, ■ ■ ■ , tn) is entire, 

(3.13) g„(--- ,t, + ^,---) =^lQn(--- ,t„---)al {l<j<n), 

(3.14) Q„(. . . , t^. + I, . . . ) = a|g„(- ■ • , t„ ■ ■ ■ (1 < J < n), 

(3.15) Rjj+i(tjj+i)Qn{- ■ ■ , tj, ^jr+i) ■ ■ ■ ) 

= Qn(- ■ ■ , tj+i, tj, - ■ ■ )Rjj^i{tjj+i), (2 < j < n - 1), 

(3.16) Qn(ti, ■ ■ ■ ,tn-l,tn)'Pn-i^n = Qn-2 (^1 , " " " ! ^n-2)'^ri-l,n) 

tn-l=tn+ri 

(3.17) trig„(ti,--- ,t„) = 0, 

(3.18) tr„Q„(ti, ■ ■ ■ ,tn) = Qn-l(ti, ■ ■ ■ , tn-l). 

These relations are verified in a way similar to those in j3] . The derivation of (j3.13|) - 
()3.14|1 rests on the transformation laws of the Xan\ which we discuss in Appendix 

El 

From the properties ()3.12|1 . ()3.13j) . ()3.14|1 . Qn can be written as 



(3.19) n^i(2ti,) X Q„(ti,--- ,t„) 



i=2 

3 n 



E-TT 6'Q,, + l(2t 



O 1 ■ ■ ■ U„ 



ai,---,o„=0 



with some Kq,^ ... G C. Terms with ai = are actually absent in the sum, in 
accordance with ()3.17|1 . For convenience we set ,02,0 = 0. Note that ()3.13j) . 
()3.14|1 and the translation invariance imply 

(3.20) i^a„,-,ai = unless Yl'j=i^j = (0)0)- 

By induction, assume Qm = for m < n. We are going to argue that Qn is then 
determined up to a multiplicative constant (see Lemma f3. II below) . 
By ()3.16|) . the induction hypothesis and ()3.15p . we have 

(3.21) Q^(... ,t^.,t^.^,,...)p-^.^^ =0 (2<j<n-l). 

By fl3.18|) we may also assume 
(3.22) 

f^an,-- ,ai — unless a„ 7^ 0. 
Quite generally, consider a matrix of the form 

Ui,2{u,v)= ^ Kp^- —- 7-rcr,a^. 
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Then the relations 

(3.23) Ri2{u - f )f/i,2(M, v) = t/i,2(f , u)Ri2{u - f ), 

(3.24) Ui,2{u + v,u)V^,2 = 0, 

are equivalent to the following relations for the coefficients Kba- 

l^b,a — l^a,b = i(yt^c,0 + «^0,c), 
3 

0=0 
3 

^ ^ Jaf^a,a = 0. 
a=l 

Here a, b, c are cychc permutations of 1, 2, 3, and Ja, Jbc are as in ()2.10|) . ()2.1ip . The 
above relations have the same form as those derived from the quadratic relations 
()2.8|) . ()2.9|) and from the Casimir elements ()2.13|) . respectively. Consider the quotient 
A of the Sklyanin algebra modulo the relations that the Casimir elements are zero. 
This is a graded algebra, 

oo 

= An ■ 

n=Q 

From the above considerations one easily concludes that there exist three linear 
functionals (a = 1, 2, 3) on An-i such that 

which satisfy the additional condition 

(SoA) = 0. 

In the Sklyanin algebra with generic parameter t], any monomial ■ ■ ■ Sa^^ can 
be reduced to a linear combination of ordered monomials Sq° S^^ S'^'^ 82^ with t^i, t'2 G 
{0, 1}, by using the quadratic relations and Casimir elements (PBW basis) [TUl lllj . 
Together with 1)3.201) this means that each functional Ka is defined by one constant, 
that is, 

(3.25) Ki{S;'~^Si), K2{S;''^S2), KsiS;'-^) for n even, 

K,{S^~''S2), K2iS^-''Si), K-i{S^-^SiS2) fornodd. 

There remain three coefficients. In order to finish the proof of the Cancellation 
Identity, it remains to show that these coefficients vanish. 
In addition to (|3.15|) . we have also the relation 

Rl2{Xl2)Qn{h, t2, ■ ■ ■ )^12(Al2)"' = Q^'^h, ^1, " " " )• 

(2) 

The poles of the R matrix in the left hand side are not the poles of Qn ■ This entails 
the relation 

Pr2Qn(tl,t2,---)n2(l) =0, 
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which can be rewritten in terms of functionals Kjq^ clS follows: 

(3.26) 2K,{ASo)-iil + Jbc)f^c{ASh) + iil-Jbc)f^biAS,) = VA G ^-2- 

These equations can be viewed as a system of linear equations for three constants 
()3.25|) . Certainly, these equations are not explicit since for every A we have to 
perform the procedure of reducing to PBW form. This huge system of homogeneous 
linear equations does not allow us to prove that the constants in question vanish; 
rather they reduce them to one constant. Let us explain this point. First, it is clear 
that the equations (|3.26p correspond to the following relation in An 

2SoSa - i{l + Jbc)SbSc + i{l - Jbc)ScSb = 0, 

obtained by solving ()2.8|) . ()2.9|) for SoSa- So, it is easy to see that all our equa- 
tions including ()3.26p are satisfied by the following construction. Consider a linear 
functional k on An such that k {Sa^ ■ ■ ■ Sa„) = unless J2]=i — (0' 0)' 

k{ASo) = 0. 

Then all the requirements are satisfied by 

Ka{A) = K{ASa). 

On the other hand the number of solutions to the system of linear equations ()3.26|) 
for three constants cannot be bigger for arbitrary rj than it is for rj = 0. In the latter 
case the algebra is commutative (see Appendix [XJ notably ()A.4|) ). and the equations 
(j3.26|) become 

k':^ (SbA) = Kf (SaA) 

with additional condition n'^ [SqA) = 0. It is easy to see that this gives a system of 
three equations for three constants whose rank equals 2. 
Thus we come to the conclusion: 

Lemma 3.1. Under the induction hypothesis, we have 
where k is a linear functional on An satisfying 

n 

K {Sa„ ■ ■ ■ Sa^) = unless dj = (0, 0) 

k{ASo) = k{SoA) = 
and as such is defined by one constant: 

k{S^) for n even, 

K^S^^"^ S1S2) for n odd. 

Unfortunately, we were not able to show that this remaining constant equals zero. 
The problem is still open. 
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4. Examples 

In this section we write down the Ansatz in the simplest cases n = 2,3. We also 
consider the trigonometric limit. 

4.1. The case n = 2. In the case n — 2, ' (^i)^2 

) can be readily found from 

()2.21|) . The function h2{ti,t2) is given as follows: 



1 1 ^ 

h2ih, t2) = -- Ha+l{2U2) O'^ ® 0\ 



4 4[ti2.„^, 
where 

This gives the formula for the nearest neighbor correlators of the inhomogeneous 
chain: 

(4.1) «^2") = -|||^'^+i(2^)' 

where a = 1, 2, 3 and t = ti2- Noting that Ha+i{2t) is odd in t, we obtain {<jf<J2) = 
—9i{2ri)H'^_^_i{0)/9[ in the homogeneous limit t — > 0, or more explicitly we have 

p 9"^ 

(4.2) ((JiCTs) 



W['9^{2r,) 

X (^2^:+i(0)e,+i(2r/) +0,+iel+i(2r/)|- -47rz0,+i0„+i(2r7)^^ ..^(O). 

Let us check the formula (j4.2j) against known results. As is well known, the XYZ 
Hamiltonian is obtained by differentiating the transfer matrix of the eight-vertex 
model 

Ti(t) =tr(i?oL(t)---i?oi(t) 



as 



(4.3) T,(0)-^n(0) = ( E<(0)«+i ) + Lv'M^ 

j=l \a=l J 

where L is the length of the chain, and we have set R{t) = J2a=o'^oi{t)o'" ® cr". 
As it was mentioned already, the R matrix (j2.H) is so normailized that in the ther- 
modynamic limit L — s> oo the free energy per site of the eight-vertex model is 0. 
Therefore, taking the ground state average of (j4.3p . we obtain 

3 

(4.4) ^r«a2^) = -/°, 

a=l 
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with J° = w;(0)^i(2ry)/^;. Explicitly we have 



(4.5) / =— (a =1,2,3), 

The average over the normalized ground state has the property 5{Hxyz) = {^Hxyz), 
where 6 stands for the variation of the coefficients Hence we have in addition 

^ dl\ „ 9/0 



a=l 



The nearest neighbor correlators {cr 1(^2) are completely determined by the lin- 
ear equations ()4.4p . ()4.6|) . ()4.7|) . Using Riemann's identity and the heat equation 
AmdO a{t\T) / dr = d'^9a(t\T)/dt'^ , one can verify that our formula ()4.2|) indeed gives 
the unique solution. 



4.2. The case n = 3. Let us proceed to the next case n = 3. Written in full, 
^3(^1,^2,^3) reads 

a+/3+7=0 



The coefficients /^^/^^ly are given as follows: 



O2 { ^4(2tl3)^3(2t23) ^ , , ^3(2^13)^4(2^23) 
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-'0,1,1 ~ 

'^'0^^ - ^ 1 ^4(2,)^3 ""^^''"^ + ^3(2,)^4 ""^^'^ 
^ft'o =2[tl3]Mi^2(2tl2), 

.(1,2) . r. i g4(2t23) WMM^r, m(9f ) 

r(l,2) _ , .^ /M4(2r7)^3(2ti3) . „ ^ r ^^si^^ rr ( ^ 

^'^■^ " i^^2(2r^)^^4 e,{2rt) [^^3]^4(2tl2) - [^13]^^i^2(2tl2) 
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.(1,3) _ .(1,3) _ 
-'0,1,1 ~ -'1,1,0 ~ 



''''' ~ 9Mi2vm2rj) i^i^ni23ji^2(2t 



^2 f ^4(2tl2)^4(2t23) ^ ... , , ^3(2ti2)g3(2t23) ^ . 



62(27]) { e,i2r])e, esi27])e: 

r(l,3) _ / .X / ^2 ^3(2^12) r N 1 ^4 ^3(2^23) z^, n\ 

r(l'3) ^ ■^ /r+ 1 ^3 ^4(2t23) rr /o+ A d2 9i{2ti2), ITT try. \\ 
'^'^'^ = ^-'^ l['-]^^^''^^2tl3) - ^^^Mi^3(2t,3)| , 

& = o> 

,(2,3) _ O2 j9j2h2)mt^ , , g4(2tl2)g3(2ti3) ^ . 

- ^ 1 ^^3^4(2r,) ""^^''^^^ + ^3(2^)^4 ""'^^'"'^ 

4S=2[tl2][tl3]i^2(2t23), 

r(2,3) _ / / ^3(277)^4 ^2(2ti3) 1 ^2(2tl2) . .1 

= H) {[hsf-^^H,i2t2s) - ^3g4(2r^)g2(2ti3) ^ J ^ . 



The rest are given by the cychc change 1^2^3^1of the indices 7 in 
^a'^\ with the change 2^3^ 
of the inhomogeneous chain are 



'with the change 2^3^4^2of the indices in 6a and Ha. The correlators 



With the abbreviation H'^ = H'^{0) and H'" = H'^'{0), we obtain a new formula 
for the next nearest neighbor correlators for the homogeneous chain 

,aa. ^ _ 1 gl(2r^) f 0!A^V)Ol+, + OlA^VWhl rr. 

^ ' A e[ \ 0,+i^^,+i05+i(2r7)e,+i(2r/) 

+ 1 J ^^.+i(2^)t^^.+i(2r^)^^^,+i ^^.+1 ""^^^^ 

We have in addition 

{ay,a^,) = 0, (aja^ag") = 0. 

In both formulas, (a, b, c) = (1, 2, 3), (2, 3, 1), (3, 1, 2). 

M. Lashkevich communicated to us a program for numerically calculating corre- 
lation functions from the integral formula of [23 12^1 • For n = 2 and n = 3, we found 
agreement between their results and ours to within the precision 10~^. 
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4.3. Trigonometric limit. Finally we briefly touch upon the trigonometric limit, 
and discuss how various quantities which appear in (j4.1|) are related to the trigono- 
metric counterpart. 

First we consider the limit to the massive regime. For this purpose, it is convenient 
to rewrite the R matrix in terms of the parameters t',rj', t' in ()2.30p as 

m = Y^,{u®u)r'{t'){u®ur\ 

where [t']' = 9i(2t'\r')/9i{2r]'\r'), r'{t') is obtained from ()2.2|) by replacing t.rf.T by 
t\rj',T\ and ^ = ^| hmit r' — +ioo while keeping A = t'/V and 

u = Irj flxed, the R matrix tends to 

(4.8) i?xxz(A) = pxxz(A) .,'';f,f^^ . 

[A + v\xxz 

In the above, [A]xxz = sinvrz/A/ sinvrz/, and 

I A /^^ 1 /sin (A + l/2)7rz/ r. „ 

4.9 rxxz A = - TT^^ ® ^ 

1\ sm7rz//2 

112 2 COS (A + l/2)7rz/ 

cos7rz//2 

I A 1 n\ (\\ ^ )oo(C )oo 

(4.10) Pxx.(A)=-C,^,^.,,^(^,)^, 

where C, = e^"^"^, q = e'^'"^, (x)oo = njlo(-'- ~" l'^''^)- It is easy to see that 

Stt 8^ ~ / a \ 

(^l{t> -. t^(A), UJ2{t) UJ{X), UJ3{t) 0, 

sm Tcu sm vrz/ 

where the functions uj{X),oj{X) are given by |4j, eqs.(13.2)-(13.5) for the massive 
regime. Hence the limit of ()4.1|) becomes 

q + q-' n^^^^±il^_~n^^ 



(4.11) hm 



^lirn^ (-1^2) = (^1-2) = 2 ( (^^^-(^) + {q - q~'^)iC - C-'f ^^^) 



which reproduces the formulas in the massive regime (see |4], Example in Section 
3). 

Second let us consider the limit to the massless regime. We set 

1 u uX 

T = r, 7] = r, t = r 

ni 2711 27ii 

for a constant u [0 < u < 1) and take the limit r | with u and A flxed. 

The limit of the R matrix is given by the same formula (j4.8|) - (j4.9p . with pxxz(A) 
being replaced by 

52(-A)52(l + A) 



Pxxz(A) 



52(A)52(1-A) • 
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Here 5*2(0;) = S2{x |2, l/u) signifies the double sine function. In tfie limit we have 
rcJift) : u;(A), ruj2{t) : ^^(A), 

Sm TTU Sm TTU 

where now uj{X) and lj(A) stand for the functions given by j4j, eqs.(13.2)-(13.5) for 
the massless regime. Moreover we have 

^ f d d 9 \ , 

' \dt '^Ot] J ^ 

From the formulas above, we see that in the massless limit the function /i2(^i)^2) 
tends to the solution /i2(Ai, A2) of the reduced qKZ equation given in [1]. 

Appendix A. Existence of TrA 

For every finite-dimensional representation of the Sklyanin algebra A, we can 
define the trace, which is a functional on A whose main property is cyclicity. In order 
to formulate our Anstaz for correlation functions, we need an analytic continuation of 
this functional with respect to the dimension. We denote this analytic continuation 
by TrA^, where A E A and X = k + 1 for 7i^''\A). In Section ESI we presented the 
formulas for TrA^Q, TrA^^ (a = 0, 1,2,3). In this appendix, we discuss the general 
case. In fact, we prove that for generic parameters Ji, J2, J3, the definition of TrA^ 
for general A E A can be reduced to these known cases. 

Consider the polynomial ring F = 'K[Kq,K2] with K = C(Ji, J2, J3). Here, we 
consider Kq, K2, Ji, J2, J3 as variables, whereas they are parameterized by r, r] and 
A in Section 121 and Appendix^ We use the parameterization in order to define finite 
dimensional representations. The discussion in this appendix is mainly concerned 
with the algebraic relations in the Sklyanin algebra only. 

We denote by A the Sklyanin algebra defined over the field K. It is a graded 
vector space, 

A = ©^^oA„, dimA„<oo. 

Multiplication by the central elements ()2.13|) endows A with an F-algebra struc- 
ture. Suppose we try to define some F-linear functional Tr on A which satisfies 
cyclicity Tr (AB) = Tt{BA). Then the question is, for how many independent ele- 
ments of A this functional should be defined. In other words, describe the F-module 



H = A/A' 



where 



Note that H = ©^=oHn where H„ = An/K, K = ELoi^a, A„_i]. 
We prove 

Theorem A.l. The F-module Ji is a rank 7 free module generated by the monomials 
(A.l) (mi)i<i<7 = (1, 5*0, 5*1, 5*2, 5*3, Sq, S^). 
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The F-hnear independence of these elements follows from fl2.20|) - ()2.23p . Indeed, 
suppose there is a relation Yli=i = with Cj G F. The sum over elements of even 
degree and of odd degree must vanish separately. Specialize Jj to the value (j2.1ip 
with ?7 Q + Qr, Imr] > 0, and take the trace of both sides on the representation 
^Qj. ^ g Z>o. By Lemma fC .21 it follows that Q = except for i = 3, 4, 5. To see 
that the latter vanish, it is enough to apply the automorphisms l^, (see ()D.1|) . ()D.2|) 
and two lines above) and take the trace. 

Let us prove the spanning property. 

Consider the tensor algebra T over the field K generated by four independent 
variables Sq, Si, S2 and S3. Set R = T[Kq, K2]. It is a graded algebra: R = ©^^o-^^fi' 
where we have dimxRn < 00. We have the isomorphism of K- vector spaces 

(A.2) H„ ~ R„,/ R„-2 ([5o, Sa] - iJbASbSc + S,St)) 

+ Rn-2(['S'b, Sc] — i{SoSa + SaSo)) + ^^[Sa, Rn-l] 
3 3 \ 

+R„_2 Si - Ko) + R„-2 JaSl - ^2) . 

«=0 a=l / 

The K- vector space R„ is spanned by the monomials of the form K^^ Sa^ ■ ■ ■ Sai 
where 2mo + 2m2 + Z = n. The relations which define H„ in ()A.2|) are linear relations 
for these monomials. For each n the coefficients of these linear relations form a 
matrix Ain with entries in K. 

The spanning property is clear for n = 0, 1. Suppose that n > 2. Divide the set of 
monomials of degree n into two groups: the first group is the monomials such that 
the part Sa^ " " " •S'q,, is equal to one of (1 < z < 7) and the second group is the 
rest. The matrix Mn is divided into two blocks = {Ai'n, -M'^) where (resp., 
J^n) corresponds to the first (resp., the second) group of monomials. It suffices to 
show that the rank of A^" is equal to the cardinality of the second group. 

The proof of this statement exploits the classical limit: 



(A.3) Ja = l~e'ja, 

We introduce new variables Sa (a = 0, 1, 2, 3) and fco, fci by 
(A.4) So = eso, Sa = sa (a = 1,2,3), 

(A.5) Ko = ko, K2 = ko- e^ki. 



In Appendix |Bj the classical limit is taken as r/ ^ instead of e ^ 0. In this 
appendix, we avoid the parametrization by r, t] and A in order to simplify the 
argument. 

In the limit e — > 0, we have 



(A.6) [sa, Sfs] = ie{sa, S/3} + 0{e^), 

where the Poisson bracket is defined by 

(A.7) {So, Sa} = 2jb,cSbSc, 

(A.8) {Sb,Sc} = 2SoSa. 



SKLYANIN ALGEBRA 21 

Here ja,b = —jb,a = ja — jb and (a, b, c) runs over cyclic permutations of (1, 2, 3). In 
the classical limit, the variables Sa become commutative. Let K'^' = C(ji, 72,^3) de- 
note the field of rational functions in ja , and let F'^' = K^^ [ko ,ki], A'^' = K*^' [sq , Si , S2 , S3] 
denote the polynomial ring in indeterminates fco, ki and Sq, (0 < a < 3), respectively. 
The Casimir relations ()2.13|) become the algebraic relations 

(A.9) sj + sl + sl = ko, 

(A. 10) si + jisl + j2sl + j^sl = ki 

in K'=^[so,Si,S2,S3, /co, A;i]: 
A^^ ~ K'^^iso, si, S2, S3, fco, ki]l{sl + S2 + S3 - fco. So + jis\ + J2S2 + J3S3 - ki). 

This makes A'^^ an F'^'-algebra. The algebra A'^' is graded as well: A'^^ = ©^g^n- 
Let (A^^)^ C A^' be the limit (in the appropriate Grassmannian such that we con- 
sider Sa as commutative variables) of A^ as e — > 0. Then, we see X]a=o{'^a' ^n-i\ <^ 
(A'^^)^ from ()A.6j) . In order to show the spanning property of (mj)i<j<7 in H, it is 
therefore sufficient to show the spanning property for 

(A. 11) (mf )i<i<7 = (1, So, Si, S2, S3, Sq, S3) 



in H'^^ = ©^=oHn\ where 



(A.12) H^> = A:;V^K,AtJ. 

a=0 

In conclusion. Theorem I A. II follows from 

Proposition A. 2. The F'^^-module H"^^ is a rank 7 free module generated by the 
monomials (jyl.iijl . 

Proof. The F'^^-linear independence of mf follows from the same argument as in the 
quantum case. In place of trace on V*-'^^ we use the non-degenerate pairing between 
cycles and co cycles given in Appendix iBl 
We prove the spanning property. Define 

^ d . d . d 

Vo = J2,3S2S3^ hj3,lS3Sl^ hjl,2SlS2^, 

osi as 2 c'S3 

Odd 

Va = -Jb,cSbScT: h SqScT^ SoSfo — , 

OSq OSb OSc 

where (a, 6, c) = (1,2,3), (2,3,1), (3,1,2). We have V„P = |{s„,P}. These are 
F'^^-linear. 

We want to show that modulo Yl^a=Q'^ a-^^^ any monomial s^° s^^ s^'^ s^^ can be 
reduced to an element in F^^ ■ 1 + YL=q ■ + ELo ■ s^. Set 

0<a</3<3 

3 



0=0 
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We also denote H'^^™ = and A^^™ = A='. 
Since 

3 oo 
Q=0 n=0 

by using ()A.9|) and ()A.10|) the above statement follows from the following: 

3 3 

(A.13) H=i'-^i = F'%so + J2^'\h-jako)sa + J2^'\h-jako)sl. 

a=l a=l 

Note that 

7 't. — 2,- 2,- 2 

"^1 JarSO — '5g + j6,a'S^ i Jc,aScy 

and therefore we have fco^o; (^i ~ jako)sa, {ki — jako)sl G A^^'~^'. 

Let us prove ()A.13|) . Suppose that a monomial m = Sq''s"^S2^S3^ G H^^'"^' is such 
that tl{o|na G 2Z + 1} > 2: e.g., no,ni G 2Z + 1. By using the Casimir relations 
()A.9|) and ()A.10|) we can replace Sq and with S2 and S3. Therefore, we have 

m G SoSiF'='[s2,S3]. 

Since 

^2{S2S3) = kSoSiS^S^ ^, 

we have 

m = 0. 

Next consider the case tl{a|na G 2Z + 1} = 1. Suppose no G 2Z + 1. Then, we 
have 

mesoslF^'[slsl] + SoslF^'[sl,sl]. 

Since 

Vi(si44) = /sosi^-'^""' - ksos\4^\s'^-\ 

the monomial m belongs to K.'^^sqsI-' C H'^''^^' where 2j + 1 = no + ni + 77-2 + na. 
Similarly, we see that 

kiso G K''soS2^ 

Since /cq'^o is a non-zero element in H'^''"^', we have 

m G K^^kisQ. 

The case n^ G 2Z + 1 (a = 1, 2, 3) is similar. 

The remaining case is no,ni,n2,n3 G 2Z and tl{a|^a > 0} > 2. We have degm > 
4. Note that 

0<j<i<3 
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We have the following relations in 



VI \ ■ 22|22 22 

V2(soS3Si) = -jg^iSgSi + SgSi - s\s\, 

V3(soSiS2) = -h,2s\s\ + s%s\ - s\s\. 

Therefore, we have dimK;ciH'^^4~^' < 3. On the other hand {k\ — jako)Sa £ for 
a = 1, 2, 3, and they are K^^-linearly independent. Therefore, we have 

Observe that for (a, b, c) = (1, 2, 3), (2, 3, 1), (3, 1, 2) we have 

-Ln , e"0+2 ria n^-l ric+l _ „no+2 ria n^+l ric-l 
-t-/ifeiQ b^ ILcOq b^ b^ 

We can increase the power in sq by rewriting Sq°s"''sJ^'''''^s"'=+^ in terms of Sq°~^'^ s]^"^ s^''~^ t 
and So°+^s]^-sl'^^s';^^~\ Thus, we see that the K^^ vector space H^^^^-" (j > 3) is 
spanned by 

sl'sf'-'^ (a=l,2,3;l<rf<j-l). 
On the other hand, the space H'^'2j^' contains K'^'-linearly independent elements 

(A.14) k^.kf'^h - jah)sl (a= l,2,3;0<rf<j-2). 

Therefore, the elements ()A.14|) span 

In the rest of this appendix, we explain a mathematical background of Proposition 
IA.21 which is the de Rham cohomology of the affine algebraic variety defined by the 
two quadrics ()A.9|) and ()A.10|) . Although our proof is independent, the statement 
of Proposition IA.2I is closely related to a result of K. Saito. 

We set 

M = {(ji, J2, J3) e C'\ja + Jb for a ^ 6}. 

Let 

(A.15) y?:X = C^xM-^r = C2xM 

be the mapping such that if> = {ipi, ...,(p5) and for (s, j) = (sq, Si, S2, S3, ji, j2, is) e 
X M we have 

(A.16) (/?i(s,j) = sl + sl + sl 

(A.17) ¥'2(s,j) = si + jisl + j2sl + jssl, 



(A.18) (/?3(s,j) = Ji, 

(A.19) (/?4(s,j) = J2, 

(A.20) (/^5(s,j) = J3. 

The critical set C C x M of this mapping is given by the equation 

dip\ A ■ ■ ■ A dip5 = 0. 
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We have 

3 

C= |J{(s,j)|s^ = 0for/3 7^a}. 

We have the commutative diagram, 

X ^ C 

Y D D, 
where the discriminant set D is given by 

(A.21) D = {{ko,h,j,,j2,j3) eC^ X M\A{ko,h,j,,j2,j3) = 0}, 

3 

(A.22) A = kollih-jako). 

a=l 

The inverse image ip~^{0) is called the simple elliptic singularity of type [221 ■ If 
y does not belong to D, the inverse image Xy = ip~^{y) is a non-singular affine 
complex surface and is called a smoothing of the singularity. The mapping 

(A.23) ^\x-^~HD) : X - ip-\D) ^Y-D 

is a locally topologically trivial fiber space, and the homology group is of rank 7: 
H2{Xy,7j) = l7 . In Appendix iB} we construct cycles in H2{Xy,'L). (In jSI], K. 
Saito defined the extended affine root systems. The homolog y group H2{Xy,'Zi^ is 
isomorphic to D^'^^ in his classification. We have not identified our cycles in D'^^'^K) 
Let Q^x t>e the sheaf of Ox modules consisting of germs of holomorphic p forms 
on X, and ^^x/y quotient sheaf 

5 

(A.24) Q^x/Y = ^x/ ^ 

i=l 

The relative de Rham complex {Q*x/y^ dx/Y) is defined by the commutative diagram 



dx/Y 








T 


T 


d 


n^x - 





This is an exact sequence of Ox modules. The following is K. Saito's result |22]- 
Theorem A. 3. The cohomology group 

H\^^{Q'x/y)) = Ker(^.(^]^/^) '^4" ^.(^]^/^))/Im '^4" V*i^'x/y)) 

is an Oy locally free module of rank 7. 

We connect the above algebro-geometric setting to ours. The following proposition 
is a corollary to Theorem IA.31 Here we give a proof in the line of this appendix 
without using Saito's result. 
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Proposition A. 4. Consider a complex of F"^^ -modules: 

= A^irfs,, A--- Ads,^, 

0<Qi<-<ap<3 

i=i,2 

The cohomology group H^^'~^' =^ Ker(Z^ —>■ Z^) /lm{Z^ —> Z'^) is a F*^'-/ree module 
of rank 7, where the action of {resp., ki) is given by the multiplication of ipi 
[resp., ip2). 

Proof. The key idea of the proof is to identify A^^'~^' with Z^. 
Let us introduce a holomorphic section u of the sheaf ^x/y 



x-c 



dsi A ds2 ds2 A ds^ ds^ A dsi 

S0S3 SqSi S0S2 

dsQ A dsi dso A ds2 dso A ds^ 

j2,3S2S3 js.lSsSl jl,2SlS2 

Consider the F'^'-module 



^cl def 



,j2,h] ^PO, Si, S2, S3,Ji,j2,j3\UJ 



We have a canonical isomorphism of F -modules 

sending P G A'^^ to Puj G A'^^ It is easy to see that 

A^"-^'^ = Z^. 

We have 

(A.25) d{Pdsa) = -Va(P)c^ 

for P G A'^'. From ()A.25|) we have 



Im(Zi^Z2)= IJ^V.A^' 

v«=0 



We have already constructed the F^^-bases of the modules A'^'/ Yl^a=o ^a^^^ ^^"^ 
A'^i'-^V ELo '^"^'=1. We will construct a basis of Ker(Z2 ^ Z3)/Im(Zi ^ Z^). 
First, observe that 



3 00 



= K'^^s^rfsi A ds2 A rfss + 5^ K^^s;^rfso A dsf, A ds^ 



n=0 a=l n=0 
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where (a, 6, c) = (1, 2, 3), (2, 3, 1), (3, 1, 2). Prom this we see that 

0<a</3<7<3 m,n>i 0<a<l3<3 

">i \l<a^6<3 l<a<2 / 

l<o<6<3 

We know that the F'^'-module Z'^/lm{Z^ -i- Z^) has the free generators kq = koSoLU, 
i^a = (^1 - jako)saUJ (o = 1,2,3) and p„ = {ki - jako)sluj (a = 1,2,3). For each 
degree e > 3 and color c = 0, 1, 2, 3, we want to construct elements of degree e and 

color c in Ker(Z^ Z^) by taking K'^^-linear combinations of k^k^K^ {a = 0, 1, 2, 3) 
and kl^k^pa (a = 1, 2, 3). We set po = ^oSq^- We have the relation Po = Pi + P2 + Ps- 
A straightforward calculation shows that for e = 3, 4 we have none; for e = 5 we 
find ^0 = ^o'^o for color and — {ki — jako)Ka for color a — 1,2,3; for e = 6 we 
find two color elements: 

^4 = (^1 - jiko)pi - {ki - j2ko)p2 - ji,2koPo, 

= (^1 - j2ko)p2 - {ki - j3^o)P3 - J2,3^0P0- 

For e = 5, 6 the above elements span the degree e cohomology classes. For e = 8, 
we have 4 obvious elements ki^j {i — 0,1; j — A, 5), and in addition, we find 

^3 3 
^6 ^ 4^0 E E + ^0 E^^i ~ jako)pa- 

a=l b^a a=l 

These 5 elements span the degree 8 cohomology classes. 

Finally, we show that Ca (0 < a < 6) are the generators of the F^'-module 

Ker(Z2 4 Z^)llm{Z^ 4 Z^). For odd e = 2n + 3 > 7, let us consider the 
color case. The cases of other colors are similar. The degree 2n + 3 and color 

space in Z^/Im(Z^ — > Z'^) has the K'^^-basis consisting of n + 1 elements kik^~^^Ko 
(0 < j < n). We have 

d(A;i Ko) = Ssl^'dsi A ds2 A ds3 G Z^. 

Therefore, the K^'-dimension of the degree 2n + 3 and color subspace of Ker(Z^ 

Z^)/lm{Z' A ^2) is n. Since the K^'-linearly independent elements k^ki ■'^0 
(0 < J < n — 1) belongs to this subspace, they span the subspace. 

For even e = 2n + 4 > 10, the degree 2n + 4 subspace of Z^/lm{Z^ Z^) has the 
K^^-basis consisting of 3(n + l) elements k^k^ ''^ pa {0 < j < n;a = 1, 2, 3). A simple 
calculation shows that the elements d{kiPo), dlk^pa) {a = 1,2,3) are K'^'-linearly 
independent in Z^. Therefore, the K'^'-dimension of the degree 2n + 4 subspace 
of Ker(Z2 4 Z^)/lm{Z^ 4 Z^) is 3n - 1. On the other hand, we obtain 3n - 1 
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independent elements of the subspace from the degree 6 elements P4,p5 and the 
degree 8 element pg. We conclude that they span the subspace. □ 

We have explicitly constructed free bases for the spaces H'^^'"' D H'^''"^' D H^^'"^'. 
In 122], similar spaces are studied for hypersurface singularities E^, Ey, Es- The 
results in this appendix supplement some part of his construction in D^, where the 
elliptic singularity is given by a complete intersection of two quadrics. 

Appendix B. Cycles and integrals in the classical limit 

In the classical limit, the defining relations of the Sklyanin algebra turns into 
()A.9|1 . ()A.10|) which define an afiine algebraic surface S C C^. In this section, we 
study the classical limit of the functional TrA and identify it with an integral over a 
cycle on S. 

Sklyanin's formulas for the representation ()2.16|) gives an explicit uniformization 
of S. Consider the limit 

T] ^ 0, A — s> oo, ri\ = fi finite. 

In the right hand side of (j2.16p . we replace rjdu by — 27rif where v is the variable 
canonically conjugate to u, 

2it{v, u} = 1. 

In the limit, the formulas for sq = So, Sa = rjSa (a = 1, 2, 3) tend to 

, , ^„+i(2M-p)e-2™-^,+i(-2M-p)e2™ 
(B.l) s^iu,v) = c. , 

where {u, v) G and 

1 ga+l(0) r ... 

Co = -, Ca = for a = 1, 2, 3. 

Eq. ()B.1|) provides a parametrization of the surface ()A.9|) . ()A.10|) with 
(B.2) h=('4f)\ fe=f^V#,ogft(,). 



9[ J ' \ 9[ 

For definiteness, we consider the case r G iM>o, < p < 1/2. The functions Sq 
have common periods 

ei = (l,0), 62 = (0,1), e3 = (r,2p), 

so {u,v) should be regarded as variables on C^/(X]Li ^^«)- shall consider the 
fundamental domain: 

< Re(u) < 1, < Im(M) < -r, < Re(f ) < 1. 

There are also pole divisors of Sq, at m = pj, where 

, , 1 1 + r r 

(B.3) Po = 0, pi = -, p2 = — P3 = -. 

If we neglect these divisors, there are, obviously, three non-trivial 2-cycles which 

are tori with generators (61,63), (61,62), (62,63). We denote them by 70, 71, 72. 
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One can choose the tori 71, 72 so that they do not intersect with the pole divisors. 
As for 7o, the first impression is that it hits the divisors. But actually one has to 
be very careful at this point. From the formulae ()B.1|) it follows that there are no 
singularities at ()B.3|) if u = Po,Pi and f = 0, 1/2, or m = P2,P3 and v = fi, fi + 1/2. 



Hence the actual divisors are V = V U V' with 



(B.4) 



V = Ui=o,i{fe, v)\Q< Re{v) < 1, 7^ 0, 1/2}, 
'D" = U,=2,3{fe, I < Re(t;) <l,v^^,^i+ 1/2}. 



This means, first of all, that we can modify 70 at f = and v = fi into a well-defined 
cycle without intersection with T>, as depicted in fig. 1 below. 



Re(v) 



® 



fig- 1 

Now we have another possibility. We can draw spheres 5q , 5i which have as south 
(resp. north) poles the points (po^O), (pi,0) (resp. (po,l/2), (pi,l/2)). In the 
vicinity of these points they are parallel to the w-plane, and every section of it by 
the plane Re(f) = afor0<a<l/2isa cycle around po,Pi in the w-plane. 





Re(v) 




N \ \ 







fig- 2 

These spheres do not intersect V. Similarly, we construct spheres 82, 5^ which have as 
south (resp. north) poles the points {p2,fJ>), {ps, fJ^) (resp. (j92, /^+ 1/2), (^3, /i + 1/2)). 

The homology group H2{S, Z) ~ 77 is known. We will see below that our cycles 
^0, ^1, (^2, ^3, 7o, 7i, 72 G H2{S, Z) give a linearly independent basis in H2{S, C). 
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Using our seven monomials we construct 2-forms 



1 dsi A dS2 2 2 

, Uo = Sq u, Us = Sg u, 



c"a = Sq, a = 0, 1, 2, 3. 

By direct computation we check that in Sklyanin's parameterization ()B.1|1 

(B.5) LU = duA dv. 

This formula allows us to calculate the integrals explicitly. 
Let us start with the cycles 6a- We obtain 

^/3+l(-/^) 



Sa 5a 

where eag are elements of the matrix: 



U = I Ug = 0, I Crp = 2tapCg- , 



Note, in particular, that 
(B.6) 



A 1 1 1\ 

1 1-1-1 

1-1 1-1 

\l -1 -1 1/ 



a/3 = 



for /3 = 1, 2, 3. 

For the integrals over 7^ {k = 1,2), we find 



(Ja = 0, UJ = 1, UJ = T 



Ik 



71 



72 
2 



71 



— - — ) d logea+M, 



_2,'qA+iM| (r9Mog^^.+i(/.) + 2vr.), 



^1 



72 



For 7o, we introduce 



Then, we have 



7o = 7o 



1 ^ 



/3=0 



= 

70 70 



^. = 0, /^. = 2(^^^^^^) 91og^^+,(/z), 



70 
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Using the above formulae we can calculate the determinant of the period-matrix as 
det {V) = Const • 0i{2fi) (^i(/i)^o(/i))' log ' ^'^^^ 



whence we conclude that for generic /i the period matrix is non-degenerate. 

Comparison with the quantum formulae ()2.20|) . ()2.21|) shows the exceptional role 
of the cycle 70 as r/ — > 0, 

TrA(5'„) ~ ^ / "^^ 



70 



1 /• 2 Jl (« = 0), 



(« = 1,2,3). 



70 

Actually, classical limits of all the elements of F in ()2.22j) - ()2.24j) can be found among 
integrals over the cycles 70,71,72- This fact played an important heuristic role in 
the calculation of traces. 



Appendix C. Technical Lemmas 

In this Appendix, we collect some technical matters related to Tr^. 
The first is the color conservation for Tr^. 

Lemma C.l. If A e A^"''"^ then 

Tt\A = unless {m,n) = (0,0). 

Proof. Set 

Mf){u) := f{u + ^)e'-^'\ Mf)in) := f{n + ^). 

One verifies easily that 'fi,<f2 ^ End(V'''^^) and 

y.iovr('=)(5,)o^-i = (-l)^^7r('=)(5„), 

where d = (ai, 02) G Z2 x Z2. It follows that 

trv.w7r«(A) = (-l)-trv^(„7r«(A) = (-l)"tr^(„7r('=)(A), 

whence the lemma. □ 

The next Lemma is concerned about the uniqueness of the representation of TrA. 

Lemma C.2. Assume Im?7,Imr > 0, r] ^ Q + Qr. Let gi {i = 1,2,3) be elliptic 
functions with periods 1 and t, and set f{u) = ({u)gi{u) + ug2{u) + g^iu), where 
({u) = —{l/2ni)9[{u)/9i{u). If there exists an N > such that fikrj) = holds for 
all integers k > N, then we have giiu) =0 (i = 1, 2, 3). 
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Proof. We divide into three cases, {i)gi{u) = g2{u) = 0, {n)gi{u) ^ and g2{u) = 0, 
(iii) g2{u) ^ 0. 

Since K = {krj \ k > N} has accumulation points, the assertion is evident in case 
(i). Let us show that (ii), (iii) lead to contradictions. 

In case (ii), considering f{u)/gi{u) we may assume giiu) = 1. Then fiu + r) = 
f{u) + l. Choose a point uq G C\{K U L) which is not a pole of f{u). One can find 
a sequence of integers ki < k2 < ■ ■ ■ , kn ^ oo, such that knTj tends to uq in C/L. 
Since Im?7 > 0, if we write kuf] = + bnT {an, 6„ G M), then the integer part of 6„ 
diverges. Therefore f{knri) diverges as n ^ oo, which contradicts to the assumption 
fiknV) = 0. 

In case (iii), we may assume g2{u) = 1. Set F{u) = f{u + rf) — f{u). We have 
F(kr]) =0 {k > N), F{u + 1) = F{u), and F{u + r) - F{u) = gi{u + rj) - 9i{u)- 
Hence F{u) = G{u) + C,{u){gi{u + rf) — giiu)) with some elliptic function G{u). 
From cases (i), (ii) we conclude that F{u) = 0. In particular gi is a constant. 
Considering f'{u + rj) = f'{u), we find that f{u) is a linear function. Clearly this is 
impossible. □ 

Let us sketch the derivation of the formulas for Tr^S'^, Tr^S*^. We have the 
standard functional relation 

t^'^ {t - ^vy^^ {t + lv)=Ht + V)t^'^'\t) + mt^^'^'Ht + ri) 
for the transfer matrices 

t^'\t) := trv(., (ri^^)(t - t^v) ■ --rl^fit - h)) G End(\/«^), 

where r('='i)(t) := (vrC^) O id)L(t) and 4>{t) = Uf=i[t - tj - ik/2)r]]. Choosing = 
l,ti =0 and applying ()C.1|) . we easily find (j2.2(jp . The same method is applicable 
for ()2.2H1 . We find it slightly simpler to use the difference equation for the matrices 

^a,2 ■ 

Appendix D. Transformation properties of Xn'^^ 

Let us study the transformation properties of X,!*'"'^ with respect to the shift of 
variables by half periods. For that purpose we exploit the order 4 automorphisms 
, of the Sklyanin algebra A given by 

Hri) 0i{v) Ooiv) HV) 

In terms of the L-operator they can be written as 

(D.l) i'{L{t)) = L{t + \)a\ 

(D.2) 6^ (L(t)) = L(t+^)a'^x (_i)e-^(2W-/4). 
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Lemma D.l. Let A G An be an element of the Sklyanin algebra of even degree n, 
and let 6i{t)^'^1it/r]^ = 9A,i(t) — it/v)9A,2(t) where gA,i, gA,2 o'^e as in (j2.18|) . Then 



e,{t)-^TT,^ri\A) = gA,i{t) -(- + ^) 9A, 



■it), 



ri 2ri 

where gA,?,if) = QA^iit + t) — QAAif) is an elliptic function. 

Proof. In view of Theorem lA.ll it is enough to consider elements A of the form 
m ■ 5*^, where m is a polynomial in Ko,K2 of degree {n — 2)/2. For n = 2, the 
assertion can be verified from the explicit formula ()2.21|) . 

Let It denote the two-sided ideal of A generated by Kq — A6i{tY / 6i{2viY , K2 — 
4:9i{t + r])9i{t - r])/9i{2r]y, and let zut : A ^ A/h be the projection. From ^IJ^ 
and ()2.17j) we have 



Along with (IdTTI) and (lDl2ll it follows that for i = 0,2 



Wt+r/2(^L^{K^)^ = wt{Ki) X (-l)e" 



■7ri(T/2+2t) 



The Lemma follows from these relations. □ 
Proposition D.2. The Xa^n obey the following transformation laws. 
(D.3) 



(D.4) 



2' 

In the last line the upper (resp. lower) sign in chosen for k = i (resp. k = j). In 
particular, X2*^'', Xg*^"* are elliptic functions ofti, ■ ■ ■ ,tn with periods 1, r. 
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Proof. It is enough to prove the case = (1,2). If k 1,2, this is a simple 

consequence of the transformation law of the L-operator 

Consider the case k = 1. Using the automorphism we have 



+ 2' ■ ■ ■) X [^12 + 1/2] + l/2][hp] 

p=3 




Applying Lemma iD.ll we obtain ()D.3|) with k = 1. Eq. ()D.4|) is shown similarly. 
The case k = 2 can be obtained by using the translation invariance. □ 
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